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We study the dynamics of wave propagation in nominally diffusive samples by solving the Bethe- 
Salpeter equation with recurrent scattering included in a frequency- dependent vertex function, which 
renormalizes the mean free path of the system. We calculate the renormalized time-dependent 
diffusion coefficient, D(t), following pulsed excitation of the system. For cylindrical samples with 
reflecting side walls and open ends, we observe a crossover in dynamics in the transformation from 
a quasi-lD to a slab geometry implemented by varying the ratio of the radius, R, to the length, L. 
Immediately after the peak of the transmitted pulse, D(t) falls linearly with a nonuniversal slope 
that approaches an asymptotic value for R/L S> 1. The value of D(t) extrapolated to t = 0, depends 
only upon the dimensionless conductance g for R/L <C 1 and upon klo and L for R/L 3> 1, where 
k is the wave vector and Iq is the bare mean free path. 

PACS numbers: 

I. INTRODUCTION 

Weak localization (WL) of electronic and classical waves arises from the interference of counterpropagating partial 
waves in closed loops. Its impact upon electronic conductance has been widely studied using steady-state methods 
such as magnetoresistance |2j , and can be directly visualized as enhanced retroreflection of light from random samples 
Partial waves with trajectories over a wide range of lengths contribute to these phenomena. Most of the studies 
in the past were focused on the static transport properties. For instance, for a bulk system, WL makes all states 
localized in one and two dimensions. In three dimensions, it renormalizes the diffusion constant and can lead to the 
Anderson localization transition when the Ioffe-Regel criterion k£ < 1 is met, where k is the wave vector and £ is the 
mean free path . It is of great interest to investigate the variation of WL in the time domain 0, S • Since 
the impact of WL increases with pathlength, which is proportional to time, it can dictate the dynamical behavior 
of wave propagation in a finite-sized sample. Though it has not been practical to isolate paths of specific lengths in 
studies of electronic conductance, this can be accomplished for classical waves in time-resolved measurements of pulsed 
transmission through random media 0, 0, 0, 0, 0, 0| . Time-resolved transmission measurements ^| have 
generally been consistent with diffusion theory, which predicts a simple exponential decay of the average transmission. 
The asymptotic decay rate due to leakage from the sample is 1/tjj = tt 2 D /(L + 2z ) 2 , where tjj is the diffusion time, 
D is the diffusion coefficient, L is the sample thickness, and z is the extrapolation length. Recently, however, 
Chabanov el al. -9| observed nonexponential decay of pulsed microwave transmission in a quasi-lD sample which was 
characterized via a "time-dependent diffusion coefficient," D(t). In fact, nonexponential behavior has been predicted 
for electronic systems. It was predicted that, when time is much greater than the Heisenberg time, tjj = i gtDi 
where g is the dimensionless conductance, the decay of the electron survival probability follows a log-normal behavior 
due to the presence of prelocalized states arising from some rare configurations of disorder in the medium 0, ^| . 
Recently, in a supermatrix model calculation by Mirlin ; 8j , an analytic expression has been found for the tail of the 
electron survival probability when t <c 1r- For the case of the orthogonal ensemble, it can be expressed in term 
of a "time-dependent diffusion coefficient," D(t)/Do = l — 2t/tji. Although such a linear decay of D(t) was found 
in microwave experiments 0; nevertheless, D(t) did not extrapolate to the bare diffusion coefficient at t = and 
the scaling behavior of the slope of D(t) did not agree with the above expression. Nonexponential decay of pulsed 
transmission has also been reported recently in numerical simulations in 2D 0] and in a self-consistent diffusion 
theory, which includes recurrent scattering |l7| . 

In this work, we solve the Bethe-Salpeter equation with recurrent scattering included in a manner that satisfies 
the Ward Identity 0, 0] to obtain the average time-dependent intensity transmitted through a random sample 
following pulsed excitation. Recurrent scattering is treated in the framework of self-consistent localization theory 
|l8l Hflj l and is included in a frequency-dependent vertex function, which renormalizes the mean free path of the 
system. We calculate the renormalized time-dependent diffusion coefficient, D(t), following pulsed excitation. The 
sample is cylindrical with reflecting side walls and open ends. By changing the ratio of the longitudinal dimension L 
and the radius R, a continuous transition can be made between two key experimental geometries: a quasi-lD geometry 
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with L > J?, which is commonly employed in microwave experiments, and a slab with L <C R, which is the typical 
optical geometry. For the quasi-lD geometry, we find D(t)/D — A — (2B/ir 2 gtD)t for t <C tu- The constant part, 
A, is universal, depending only on g, while B is nonuniversal and depends upon L/£ , where In is the bare mean free 
path, as well as upon g. In the limit j > 1 and L/£ ^> 1, our results coincide with supersymmetry calculations in 
Ref. [7]. For moderate values of g, g > 5, our results are in agreement with experiment |9j. For the slab geometry, D 
approaches a nearly constant renormalized value, being equal to Dq(1 — 1.03/(fc£o) 2 ) in the limit L/£ ^> 1. This is 
close to the result of WL theory for a bulk system. 



II. THEORY 



We consider a plane wave pulse normally incident on the front surface of a random sample at z = 0. We assume 
that there is no gain or absorption in the medium and that scattering is isotropic. The ensemble-averaged intensity 
{I(t,r)) within the sample is obtained from the Fourier Transform of the frequency correlation function of the scalar 
field C n (u;r) = (0 n+ (r) <^_(r)), i.e., 

(I(t,r)) = —J duexp{-iu>t)Cn(u;r), (1) 

where fl^ = Q±ui/2, f2 is the central frequency, lj is the modulation frequency and ^n(r) is the wave field at position 
r inside the sample. In fact, it was Prof. George Papanicolaou who pointed out the importance of this frequency 
correlation function in the study of wave dynamics in randomly layered media |20ll2lll22| . Cq(w; r) = (</>n+ (r) (r)), 
is obtained from the space-frequency correlation function Cq(oj; r, r') = (</>q+ (r) </>q_ (r')), which satisfies the following 
Bcthe-Salpeter equation, 

Cn(w;ry) = <<^ + (r)> <<&_(r')> 

' tin dr 2 dr 3 dr 4 (G n + (r, n)) ( G n - (r', r 3 )) 
x t/o(w;ri,r 2 ;r3,r4)Cn(w;r2,r 4 ), (2) 

where (</>n(r)) is the coherent source inside the sample, and (Gsi(r,ri)} = — ^4^177^ is the ensemble-averaged 
Green's function that represents the coherent part of wave propagation from ri to r. The complex wavevector 
k = k + ^ describes the ballistic propagation inside the random media, where k — ^ is the wavevector, v is the phase 
velocity, and £ is the scattering mean free path, which is determined from the imaginary part of the self-energy of (G) . 
In the absence of weak localization, the bare mean free path, £q , is determined from the single-scattering diagram 
via £q = 1/na, where n is the density of scatterers and a is the total scattering cross section. The vertex function Uq 
represents the sum of all irreducible vertices. Here we approximate Uq as 

C/n(w;ri,r 2 ;r3,r 4 ) = 

^[1 + 5{u, k)] 5(T! - r 2 ) 5(1! - r 3 ) «(n - r 4 ). (3) 

The first term in the vertex function with a scattering strength 4% /In represents self- avoiding paths and generates all 
the ladder diagrams that give rise to wave diffusion when L ;§> £0 23] . It is worth mentioning that in the absence 
of the second term Eq. (2) is equivalent to the radiative transport equation, which has been extensively used by 
Prof. George Papanicolaou in the study of multiple scattering in random media |24|. The second term with a vertex 
strength 4ir5(u>, k)/£n represents WL contribution to the vertex function. The presence of this term renormalizes the 
bare mean free path to a frequency-dependent mean free path, i.e., £(u>, k) — £q/[1 + S(u>, k)]. For flux conservation 
to hold, the Ward Identity requires that the mean free path £ that appears in (G) should also be replaced by 

the same £(u>,k). Since the second term represents recurrent scatterings, it is obtained by summing all maximally- 
crossed diagrams due to weak localization in the tube geometry. The assumption of point-like scattering for the WL 
contribution is justified as long as the system is far from the localization threshold, so that the renormalized mean 
free path is scale independent. In a bulk, the renormalization factor 8(u>, k) can be obtained from the renormalized 
diffusion coefficient 19], which can be written as 

' ±[l + ^GV;r,r)], (4) 



D(uj,k) D l k 2 
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where D(w,k) = v£(w,k)/3, and G is the Green's function that satisfies the diffusion equation in the frequency 
domain: 

(A)V 2 +^)G(^;r,r') = -<5(r - r'), (5) 

The diagonal term of the Green's function, G(uj; r, r'), represents the return probability of waves that travel diffusively 
inside the sample. For the geometry considered here, we solve for G by using the boundary conditions for a tube of 
radius R and length L , with perfectly reflecting cylindrical wall and open ends. The result can be written as 

G = Gi+G 2 , (6) 

with 

z = z') = f Sin2[g " (z + Z " )] (7) 



and 



1 ™ c 

G 2 (w; z = z')= ~ } J sin 2 [<7„(z + z )] 

ZTTjJnL . 
u n— 1 

-iu + D (ql + a 2 /£l) 



x In 



-tco + D {ql + 1/i? 2 ) 



(8) 



where q n = nn/L is the transverse momentum, n c = aL/irlo is the upper momentum cutoff in the z-direction, 
L = L + 2zq is the effective length of the sample and zq = 0.71044) is the extrapolation length Q. The factor a/£o 
in Eq. (8) denotes the upper momentum cutoff in the transverse direction, a is chosen as 1 in our calculations unless 
otherwise specified. The G\ term arising from all the diffusive modes that are uniform in the transverse direction, 
i.e., q± = , is responsible for the continuous change of decay rate in the ensemble-averaged transmitted intensity 
(I(t)). The existence of this term is due to the reflecting boundary conditions at the cylindrical wall at which VG = 
normal to the surface. The Gi term represents the contributions of the transverse modes other than q± = 0. Here 
we have assumed that G 2 is independent of the transverse position, which is valid when R^$> £q. This term becomes 
important when 2R/L is not small. In the limit of slab geometry, i.e., 2R/L 3> 1, G\ is insignificant and G2 becomes 
the main contribution. Eqs. (7) and (8) indicate that the renormalized mean free path is z-dependent. In order to 
simplify our calculations, we take the spatial average along the z-axis and replace the factor of sin 2 [q n (z + Zq)] by 
1/2. This simplification is valid as long as L ^> £ . The self-consistent theory of WL requires us to replace D in Eq. 
(7) and (8) by D(u),k) Using Eqs. (3,4,6-8), we finally obtain 



8(u> i k) = 6i(w,k)+6 a (w,k), (9) 

with 

5i (u,k) = — ^ V , , N , (10) 



and 



2kH L ^ \-iu + D(w, k) {ql + 1/i? 2 ) 

where N = (fci?) 2 /4 is the number of transverse modes in the tube. Eq. (11) requires R > io/a, otherwise, 
S2(u> — 0) < 0. In our calculations, we solve Eqs. (10) and (11) for S(uj, k) self-consistently with D(u, k) — v£(u>, fc)/3 = 
v£o/[3(l + 5(oj, k))}. Physically, the self-consistency of £(u, k) or D(uj, k) assures the continuous renormalization due 
to recurrent scattering paths. In the static limit and when L/£q is large, it can be shown that 8x(u> — 0) = l/(3g — 1), 
where g = AN£q/3L is the dimensionless conductance. This result implies that localization transition occurs at 
g = 1/3 or N = L/A£q in quasi-lD geometry. 
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In our calculation, we assume both the excitation intensity and the scattered intensity are uniform over the transverse 
cross-section of the tube. Eq. (1) can now be written as 



Cn(uj,z) = exp 
1 



ILUZ Z 



dz' H(lu, z — z) Ch(lu, z ) , 



(12) 



where 



H(uj, z — z') — tt dp 



exp 



(^-iRiy)vV + (--^ 



p 2 + (z - z' f 



(13) 



Eq. (12) is numerically solved for Cn(ui,z). The transmitted intensity (I(t,L)) is then calculated from the Fourier 
transform of Cq(lu,L) , from which we obtain the renormalized time-dependent diffusion coefficient by using the 
relation D(t)/D = — t d d\n(I{t)) / dt, for t > t d , where t d — L 2 /tt 2 D is the diffusion time 0- 



III. RESULTS AND DISCUSSIONS 



The time evolution of wave propagation is described by the renormalized diffusion coefficient D(t)/Do, which, in 
turns, depends on the behavior of S(oj). Thus it is important to understand the general behaviors of 6i(oj) and 82(1^) ■ 
In Fig. 1, we show the functions <5i(w) and 82(10) for a typical case of k£o = 4 and L = 20£o. For 5i(ui) we choose 
R = 3.5£q (2R/L = 0.35), which corresponds to N = 50. For 82(^1) we choose R as infinity, which corresponds to 
a slab geometry. The diffusion time in this case is td = 139£q/v, which implies that the width of the frequency 
correlation function is, 8cu = 2it/td = 0.045. From Fig. 1, a shape rise in the real part of 8\(lS) (curve a) is found 
when uj < 8m. This behavior suggests that significant renormalization of D(t) will occur when t > td. In fact, it is 
this sharp rise of 8\ at small u> that leads to the linear decay of D(t) found previously However, the behavior 

of ^2(^) is very different from that of 5\(uj). The real part of 62(^0) (curve c) increases slowly with decreasing oj. 
This suggests a gradual renormalization of D(t) and will lead to a nearly constant D(t) in the slab geometry. This is 
also consistent with the exponential decay of the transmitted intensity found previously in the optical measurements 

EEi. 

In the inset of Fig. 1, we plot <5i(0), ^(O) and 5(0) as functions of 2R/L. It is clear that Si dominates when 
2R/L <C 1. As 2R/L is increased to 1, 62 becomes the dominant term. The saturation of 62 when 2R/L > 2 implies 
that the behavior of D(t) will not change much when 2R/L > 2. 



(a) quasi-lD system 



For quasi-lD systems, we first compare our theory with the microwave experiments reported in Ref. [9|- The results 
are shown in Fig. 2 for N — 68 and L/£q ~ 6.4 (sample A), 9.5 (sample B), and 19.2 (sample C), corresponding to 
£q = 9.5 cm. The best fit is found at a = 0.5, and good agreement is found for Samples A and B. For Sample C, for 
which g = 3, however, our theory predicts a smaller decay rate, reflecting the important role of localization effects, 
which are not included here. 



The linear region of the diffusion coefficient immediately after the peak can be fitted with D(t)/Do — A — Bt" , 
where t" = 2tj'K 2 gT£,. Mirlfn's analytical result of Ref. |8| predicts A = B = 1. Nevertheless, our calculations |25j 
show that A is a universal function of g and can be fitted by a single curve A = 1.00 — 0.27/ g — 0.17/g 2 . In the limit 
of large g, A approaches 1. However, the slope B depends on both g and L/£q. For each value of L/£q, B can be 
fitted with B = Bq(L) + B\(L)/g + B2(L)/g 2 with Bq(L) being a linear function of (£74)) -1 - I n the limit of large 
L/£o, Bo approaches 0.96. Thus, our results indicates that the analytical expression of [8j is valid only when both g 
and L/£q are sufficiently large. 
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FIG. 1: 8(u>) with L — 20^o and k£o = 4 are plotted as functions of uj. Here (a) is the real part and (b) is the imaginary part of 
5i(u>) when R — 3.5£o (N — 50), while (c) is the real part and (d) is the imaginary part of 82(1^) in the limit of R — > 00 (slab 
geometry). In the inset: 8(u> = 0), 8i(u> = 0) and 82(1*) = 0) for L — 20£o and k£o = 4 are plotted as a function of 2R/L. 

(b) crossover from quasi-lD to slab geometry 

It is interesting to see how the behavior of D(t) crosses over from a quasi-lD to a slab geometry as 2R/L increases. 
A plot of D(t)/Do as a function of dimcnsionless time t' = tj td at 2R/L = 0.35, 0.5, 1, 2, 4 and 00 (slab) is shown 
in Fig. 3 for the case L = 20£o and Mo = 4. When 2R/L = 0.35, the Si term dominates and D(t)/Do decreases 
rapidly between t ~ 2tn and 15tn and gradually saturates to a constant value at later times. As 2R/L increases, the 
peak of D(t)/D near t ~ 2 — itu also increases. In the meantime, the decay of D(t)/D at later times is reduced. 
This indicates a smaller WL effect when 2R/L is increased. This is consistent with the decrease of 5(0) with 2R/L 
shown in the inset of Fig. 1. When 2R/L > 2, D(t)/D approaches a constant value. This again is consistent with 
the saturation of 5(0) when 2R/L > 2 shown in the inset of Fig. 1. 

In order to give a systematic description of the crossover behavior, we present the crossover behavior of the intercept 
A obtained from the linear fit of D(t)/Do immediately after the peak. In Fig. 4, we plot A as a function of g for 
different values of 2R/L and for L/£ = 20. The data are well fit by A = Aq + A\jg + A 2 /g 2 . Thus, at large g, 
we have A = 1 — Co(2R/L)/g, where Co is an increasing function of 2R/L. From the fit of the data (dashed lines), 
we find that Cq(2R/L) is proportional to R 2 when 2R/L > 1. This R 2 dependence cancels the R 2 factor in g in 
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FIG. 2: D(i)/Da is plotted as a function of t' = t/ro for three samples studied in Ref. Q, with N = 68, L/lo = 6.4 (sample 
A), 9.5 (sample B) and 19.2 (sample C), where Do = 42. 5cm 2 /ns. The theoretical results for the corresponding samples are 
shown in solid lines. 

the denominator. Thus, in the slab limit, we can write A = 1 — c(L)/(k£ ) 2 . This behavior can also been seen from 
the inset of Fig. 4, where A is plotted as a function of 2R/L for different values of k£ for the case of L = 20£ . A 
approaches a constant value when 2R > L for each value of k£,Q. 

(c) slab geometry 

In order to find the function c(L) in A(k,L) = 1 — c(L) / (k£o) 2 for the slab, we plot A(k,L) as functions l/(k£o) 2 
in Fig. 5 for different values of L/£o. The slopes of these curves give c(L) and is shown in the inset of Fig. 5. 
Since c(L) is an increasing of L, D(t)/Do = A decreases with L. This reflects continuous renormalization of the 
diffusion coefficient with increasing L due to the presence of longer recurrent scattering path lengths in a thicker 
sample. The fitting of the data of c(L) in the inset of Fig. 4 gives c(L) = 1.02 — 3.72£q/L. Thus, we find D(t)/Do = 
1 — l.02/(k£o) 2 + 3.72£q/ L(k£o) 2 . In the limit of large L, we recover the result of weak localization theory for the 
bulk, i.e., D/Dq = 1 - 3a/7r(H ) 2 with a = 1 
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IV. CONCLUSION 



In conclusion, we have solved the Bcthe-Salpeter equation with recurrent scattering to find the impact of WL on 
the dynamics of wave propagation through a finite-sized mesoscopic sample. The crossover from quasi-lD to slab 
geometry is studied through the change in the behavior of the time-dependent diffusion coefficient D(t) for different 
aspect ratio 2R/L. In quasi-lD systems, D(t) decays rapidly immediately after the peak of the pulse is transmitted. 
In contrast, D(t) is nearly constant for a slab. These different behaviors are attributed to two different kinds of WL 
factors which appear in each case. These results are consistent with microwave || and optical |l'2lll.'{l| measurements. 
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